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ABSTRACT 

Image quality in mosaicked observations from interferometric radio telescopes is strongly dependent 
on the accuracy with which the antenna primary beam is calibrated. The next generation of radio 
telescope arrays such as the Allen Telescope Array (ATA) and the Square Kilometer Array (SKA) have 
key science goals that involve making large mosaicked observations filled with bright point sources. 
We present a new method for calibrating the shape of the telescope's mean primary beam that uses the 
multiple redundant observations of these bright sources in the mosaic. The method has an analytical 
solution for simple Gaussian beam shapes but can also be applied to more complex beam shapes 
through minimization. One major benefit of this simple, conceptually clean method is that it 
makes use of the science data for calibration purposes, thus saving telescope time and improving 
accuracy through simultaneous calibration and observation. We apply the method both to 1.43 GHz 
data taken during the ATA Twenty Centimeter Survey (ATATS) and to 3.14 GHz data taken during 
the ATA's Pi Gigahertz Sky Survey (PiGSS). We find that the beam's calculated full width at hah 
maximum (FWHM) values are consistent with the theoretical values, the values measured by several 
independent methods, and the values from the simulation we use to demonstrate the effectiveness of 
our method on data from future telescopes such as the expanded ATA and the SKA. These results are 
preliminary, and can be expanded upon by fitting more complex beam shapes. We also investigate, 
by way of a simulation, the dependence of the accuracy of the telescope's FWHM on antenna number. 
We find that the uncertainty returned by our fitting method is inversely proportional to the number 
of antennas in the array. 

Subject headings: data analysis and techniques 
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1. INTRODUCTION 

The ATA is a "Large Number of Smah Dishes" (LNSD) 
interferometric array located in Northern California at 
the Hat Creek Radio Observatory (HCRO) of the Uni- 
versity of California, Berkeley. The telescope is designed 
to be highly effective for commensal surveys of con- 
ventional radio astronomy projects and SETI (Search 
for Extra- Terrestrial Intelligence) targets at centimeter 
wavelengths, and currently consists of 42 6-meter dishes 
with continuous f r equency coverage over 0.5-10 GHz 
(jWelch et al.1120091 : IWiUiams fc BoweillMoh . The ATA 
has a wide field of view (4.7 deg^ at 1.43 GHz), which 
allows for rapid surveying of large areas of the sky. Since 
the point-source sensitivity of wide-field surveys requir- 
ing multiple pointings is proportional to ND (as opposed 
to the total collecting area ND^), where N is the number 
of dishes and D is the dish diameter, the LNSD design 
was the natural choice for the ATA (| Sargent fc WelchI 

[Tool . 

Since the ATA is a wide- field survey instrument, the 
construction of mosaics — large images created by stitch- 
ing together snapshots from individual pointings — is of 
prime importance. In order to account for the atten- 
uation of sources with distance from the center of the 
primary beam (i.e. the pointing center), and to correct 
for this effect when mosaicking the images, we must know 
the shape of the telescope's primary beam. 
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The results of not properly ch aracterizing the beam can 
be devastating. iCornwell et al.l ([1993) discuss the effects 
of poor beam characterization on image fidelity, noting 
that deviations from the primary beam model of a tele- 
scope should be small down to or below the ~7% level. 
Wright & Corder (2008) show that if the actual primary 
beam pattern deviates from the canonical pattern by just 
a few percent, the image fidelity can be reduced by more 
than two orders of magnitude. 

Our goal is to self consistently measure the FWHM 
with science data, which allows us to monitor the beam 
shape in real time and at the current observing frequency. 
While wo rk such as the hol ography mapping study de- 
scribed in iHarp et"al] ()2010D is invaluable as a thorough 
initial calibration of the beam shape of the telescope, it is 
time-consuming, and cannot be repeated regularly. Ad- 
ditionally, holography mapping is performed by training 
the telescope on strong satellite emission, which may be 
at a frequency different from that being used for obser- 
vations. The real-time calibration work we present here 
is complementary t o such detailed studies. 

iWelch fc DeBoeii ()2004[ ) predicted the gain of the ATA 
antennas to be 
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where 6 is the angular distance from the pointing cen- 
ter. Furthermore, the beam can be approximated by a 
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symmetric Gaussian of the form 
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where O1/2 is the theoretical full width at half maximum 
(FWHM) of the Gaussian beam. 

Both the beam from Eqn. [1] and the Gaussian beam 
approximation from Eqn. [5] are plotted in Figure [Tl along 
with the residual between them. The residual peaks at 
2.5% at a gain of less than 10%; the beam shapes differ 
by less than 0.5% inside the beam's half-power points. 
These small differences suggest that the Gaussian ap- 
proximation from Eqn. [2] is an acceptable alternative to 
the more complicated beam shape from Eqn. [TJ 

Scaling the predicted Q1/2 with frequency in units of 
GHz gives 
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/(GHz) 
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where the canonical value of 6o is 3.50°. 

Equation [3] gives us an expected FWHM of the mean 
primary beam of the ATA (the average of the primary 
beams of eac h dish) of 2.45° at the ATATS frequency 
of 1.43 GHz (iCroft et al.ll2010l). and 1.11° at the PiGSS 
frequency of 3.14 GHz (|Bower et aIll2010D . The actual 
shape of the beam may not match this prediction, how- 
ever. In order to mosaic the large numbers of snapshots 
successfully, we must properly characterize the beam. 

In we explain our two methods for characterizing 
the FWHM of the telescope's mean primary beam, after 
which, in Sj3l we describe the data to which we applied 
our methods. We then report our results in fJH and we 
conclude in ij5] with a discussion of a simulation that ap- 
plies our methods to simulated data of future telescopes 
such as the expanded. 350-dish ATA and the 3000-dish 
SKA (SKA-3000). 

2. DESCRIPTION OF METHODS 

As pictured in Fig. [21 pointings used in an ATA mo- 
saic overlap significantly, having been arranged to ap- 
proximate a hexagonal close-packed grid. The distances 
between pointing centers are determined by Nyquist sam- 
pling, and the convergence toward the poles of lines of 
constant RA has been accounted for. This is a typical 
pointing grid used in many all-sky surveys with the ATA 
and other telescopes. 

Sources can appear in different locations within the 
beam in several adjacent, overlapping pointings. Bright 
sources ('-^l Jy) are sometimes detectable in as many as 
six or seven adjacent pointings, depending on where in 
the primary beam they fall in each pointing. Dimmer 
sources tend only to be detected in two or three point- 
ings. We use the multiple appearances of each source 
to constrain the beam's FWHM by comparing the flux 
densities of the source in adjacent pointings. 

To find matching pairs, we compare each detected 
source with every other source, declaring a match when 
the difference between the two sources' reported posi- 
tions is less than the typical pointing error of one arc 
minute. We accumulate 4,304 matched pairs in the 
ATATS data, and 824 pairs in the PiGSS data. After 



matching sources, we attempt to constrain the beam's 
FWHM using two methods of parameter fitting. 

This method assumes that the sky is stationary on 
the time scale of observations. Transient and variable 
sources, therefore, will reduce the accuracy of the overall 
result. Transient sources — defined as those that appear 
in a single pointing — can be rejected on the absence of 
a source match in another and will then have a minimal 
effect as long as source confusion is not important. Vari- 
able sources present a more complex problem. It is likely 
that strongly variable sources can be rejected with a tech- 
nique such as median filtering. But sources with moder- 
ate variability on the timescale of the observations will 
be difficult to distinguish from distortions in the primary 
beam shape. Fortunately, the number of variable sources 
is a relatively small fraction (less than 1% at mJy sensi- 
tivity) of the total nurnber of radio source s (| Bower et al.l 
[2007t[Crofreran[20T0tlBower et alllMl . 

We must also emphasize that the true coordinates for 
the primary beam are (Az, El), while these measure- 
ments give results for celestial coordinates (RA, DEC). A 
more sophisticated model would require translating (RA, 
DEC) to (Az. El), especially if we wanted to work with 
long tracks of data. Our current method is appropriate 
for situations where the beam is circularly symmetric, 
and when working with the snapshot observations typi- 
cal of high-cadence surveys such as ATATS and PiGSS. 

2.1. Method 1: two-point Gaussian fitting 

Assuming only a radial sensitivity dependence within 
the telescope's primary beam, we can calculate the 
FWHM of a Gaussian that we fit to each pair of matched 
sources. We represent the measured flux densities in the 
following way: 



Si — S'truo G{6i) , 



(4) 



where iStrue is the absolute flux density of the source (un- 
known to us) ; Si is the observed flux density of a source 
in a given pointing; 9i is the distance of the source from 
its respective pointing center; and G{9) is the gain of the 
telescope as a function of the distance from the pointing 
center. 

If we assume the gain is a circular Gaussian with a 
given standard deviation a, then the observed flux den- 
sities are 



Si — "S'truo exp j ^ 2 
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Solving for for a given pair of sources gives 
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Our analytical method for calculating FWHM values 
using Equation ini breaks down under two circumstances: 
when S1/S2 ~ 1 (which causes us to divide by zero), 
and when 0i/02 ~ 1 (which results in a FWHM of zero). 
(These two cases are essentially the same, since when 
61 K, 62, it follows that Si w 5*2. This is assuming a 
symmetric primary beam, and assuming that the source 
has a constant flux density and hasn't been mismatched.) 
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Fig. 1. — The beam shape from I Welch &: DeBoeil II2004I ) and the Gaussian beam approximation (left), and the residual gain after taking 
the difference between the Gaussian approximation and the Welch beam (right). The plotted beams correspond to the PiGSS frequency of 
3.14 GHz. 
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Fig. 2. — An illustration of seven overlapping PiGSS pointings. 
The crosses are the pointing centers, the solid circles denote the pri- 
mary beam's half-power radius, and the dotted circles denote twice 
the half-power radius. This hexagonal close-packed arrangement of 
pointings is typical of other all-sky surveys, including ATATS. 



This is unfortunate, because the very places where Equa- 
tion[5]breaks down are where we would expect the largest 
number of pairs to be detected: that is, where the source 
appears roughly the same distance from each pointing 
center and is not attenuated too much in either point- 
ing. 

The advantages of this method are that it calculates 
a FWHM measurement for each pair of matches; for a 
purely Gaussian beam shape, the analytical solution for 
the FWHM is simple to calculate; and the method's accu- 
racy should increase as the number of observed sources 
is increased. The disadvantages, however, also become 
apparent in our discussion above. 

2.2. Method 2: minimization 

Instead of solving for a separate FWHM for each flux 
density pair, this approach allows us to correct all of the 
data at once with one best-fit FWHM. We correct the 
measured flux densities of each pair of matched sources 
to their (presumed) true values by multiplying by the in- 
verse of the gain, which is defined according to Equation 
[21 and depends on the primary beam's FWHM. To ar- 



rive at the best-fit FWHM value, we assume a FWHM 
for the beam, correct all of the flux density pairs using 
that FWHM, and calculate the reduced-x^ {xt) using 
the following formula: 
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where v is the number of degrees of freedom, N is the 
number of matched pairs, S'i,true,i and >5'2,true,j are the 
corrected flux densities, and ASf ^ and AS"! ^ are the un- 
certainties in flux density. 

An important advantage of this method is that it, un- 
like the two-point method, it uses all of the data to deter- 
mine the FWHM of the primary beam. Additionally, we 
can use this method to fit any arbitrarily complex beam 
shape, not just simple shapes like the purely Gaussian 
shape G{6) that we have assumed above. For example, 
the method could easily be extended to fit other param- 
eters such as the beam's ellipticity and angle on the sky. 

3. DESCRIPTION OF DATA 
3.1. ATATS 

We use both methods described above to determine 
a value for the FWHM for the primary beam of the 
ATA, using da ta fro m the ATA Twenty-centimeter survey 
(ATATS; [Croft et al.ll20100 . ATATS is a multi-epoch, 
690 deg^ survey at 1.43 GHz, and is designed to search for 
transients as well as to verify the capabilities of the ATA 
to form large mosaics. ATATS i s centered on the ^^10 
deg^ Bootes deep field ([Jannuzi fc Dey 1999). The data 
were reduced using a software suite developed for reduc- 
tion of ATA data, RA PID (Rapid Automa ted Processing 
and Imagi ng of Data; 



MIRIAD (iSault et al 



Keating et ;iL|[2009f ). RAPID calls 
19951) tasks using C shell scripts. 



and performs RFI excision, visibility calibration, phase 
self-calibration, and amplitude flagging. 

The images were deconvolved using the MIRIAD task 
"clean". The number of iterations was set using an "in- 
telligent clean" algorithm implemented by RAPID. The 
flux densities of point sources in the field, and the mea- 
sured RMS noise, are used to calculate the appropriate 
number of clean components to clean down to the 3(t 
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level. At this stage, RAPID checks for the presence of 
point sources in the residual image, and if any are found, 
the number of clean iterations is iteratively increased un- 
til the residual images are consistent with noise. 

The sources we used in our analysis were 5623 sources 
from a catalog created by running the MIRIAD task 
SFIND on maps of individual pointings from A TATS 
Epoch 12, taken on 2009 April 3 (jCroft et al.l2010[) . Each 
of the 253 pointings is a snapshot with one minute of in- 
tegration time. The overlapping of the individual point- 
ings, as described in [J2]and depicted in Fig. [21 allowed us 
to detect the same source in several different snapshots. 
These data were taken as part of a pilot study; ongoing 
improvements to the array and to the data analysis rou- 
tines are resulting in continuing enhancements in image 
quality and data consistency. 

3.2. PtGSS 

The Pi Gigahertz Sky Survey (PiGSS; iBower et al.l 
is a radio continuum survey at 3.14 GHz, which 
observes the North Galactic Cap and other selected 10 
deg^ fields. These data were processed using methods 
similar to those used on the ATATS data; the images 
are restored with a circular synthesized beam of 100" 
FWHM. The data used in our analysis come from the 
individual images comprising a seven-pointing mosaic of 
the Bootes deep field. These images were obtained from 
75 daily observations, each of which included between 
two and six snapshots of each pointing, providing very 
good coverage of the (u, v) plane. Maps of each pointing 
were produced with a size of 4.26° on a side, much larger 
than the primary beam FWHM of ~1.1°. A total of 672 
sources were identified in the individual pointings using 
the MIRIAD task SFIND. 

4. RESULTS 

4.1. Two-point method 

Fig. [3] shows histograms of FWHM values calculated 
using the two-point method on the ATATS and PiGSS 
data. We truncate the horizontal axis at a FWHM of 
12° for ATATS and 4° for PiGSS, but the tail of large 
FWHM values continues as a result of the breakdown 
of Equation [S] as described in H2.1[ While the distribu- 
tions of FWHM values pictured in Fig. [3] are strongly 
peaked, they are not symmetric. We therefore estimate 
the FWHM by taking the median of our results, and 
finding the extremes that encompass 68.3% of our data. 
We find that the mean FWHMs are 2.29 ± 0.01° for the 
ATATS data and 1.10±0.01° degrees for the PiGSS data. 

4.2. minimization 

If the assumption of a radial, Gaussian beam is correct, 
the values we calculate for 5*1, true and 5*2, true should be 
the same (and if our overall flux calibration is correct, 
should be equal to the true source flux density). Fig. 2] 
shows 5*1 vs. 52 for both the uncorrected and corrected 
ATATS data. Fig. [5] shows the equivalent plots for the 
PiGSS data. 

Fig. [5] shows a plot of as a function of FWHM 
for both ATATS and PiGSS. We had 4,303 and 823 de- 
grees of freedom in our one-parameter fits for ATATS and 
PiGSS, respectively. The minimal xt values for the fits 
are quite large: for the ATATS data the value is ~21, and 



for PiGSS the value is ~10. This is n ot due to the dif- 
ferenc es between the two models from iWelch fc DeBoeJ 
(|2004f) . expressed in Equations [1] and [2j after perform- 
ing a analysis that corrected the flux densities using 
the beam shape from Equation [TJ we calculated similar 
values for xl for both ATATS and PiGSS. 

The large values of xt are due to a combination of sys- 
tematic underestimation of uncertainties in flux density, 
the fact that the beam is not a perfect circular Gaussian 
(jHarp et al. 2010), and the fact that there are probably a 
few mismatched sources. This value could be reduced by 
fitting a more complex beam shape that includes beam 
ellipticity and beam angle, which neither Equation [1] nor 
Equation [2] took into account. 

After correcting the data and performing a x'^ min- 
imization as described in i i2.21 we arrive at a best-fit 
FWHM value of 2.39 ± 0.04° for the ATATS data, and 
1.10 ± 0.01° for the PiGSS data. Fig. [7] compares our 
results with the expression for expected FWHM, given 
in Equation |3l 

Clearly we are seeing a slightly narrower beam width 
than the width predicted by [31 The bias of a few percent 
is probably due to a combination of systematic errors 
including imperfect understanding of the ATA antennas, 
the inadequacy of the Gaussian beam assumption, and 
other losses that have a radial effect. 

5. SIMULATION: ATA-350 AND SKA 

Here we use a simulation to investigate the dependence 
of the accuracy of the telescope's FWHM on antenna 
number. As the number of antennas in an array in- 
creases, both the sensitivity of the array and the number 
of detectable sources increase. Both of these factors al- 
low us to measure the FWHM of the telescope's mean 
primary beam more accurately as the number of dishes 
increases. 

We assume each simulated array has a System Equiv- 
alent Flux Density, or SEED, of -6000 (the SEED of 
the ATA-42; see Equation ITU)) . a bandwidth Ai/ of 200 
MHz, and consists of identical 6.1 m antennas. We gener- 
ate simulated sources by assuming snapshot observations 
with an integration time r of 1 minute, an observing fre- 
quency of 3.14 GHz (the same as the PiGSS data), and 
a power-law di stribution of source density as a function 
of flux density (jVerschuur fc Kellermannlll988() : 



Assuming the power-law distribution of source density 
from Equation [HI we arrive at the number of sources per 
steradian theoretically observable by a telescope sensitive 
to a given RMS flux density: 

JL^^ (9) 
An AS ' ^ ' 

where N is the number of detectable sources, Afl is the 
solid angle subtended by the telescope's field of view, 
TVo = 3 X 10^ Jy-i ster-i, = 0.01 Jy {Nq and 5*0 are 
both derived from the plot of s ource count vs. flux den- 
sity at 1.4 GHz on p. 651 of ' Verschuur fc Kellermamil 
(|198S1 ). and AS" is the RMS flux density of the observa- 
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Fig. 3. — A histogram of the FWHM values calculated using the two-point method on ATATS data (left) and PiGSS data (right). 
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Fig. 4. — Si vs. S2 for identical ATATS sources, prior to (left) and after (right) correcting for the Gaussian beam attenuation using 
the best-fit (minimum x^) value for the FWHM. The appearance of clusters of points vfith the same flux density in either the vertical or 
horizontal direction results from the fact that some sources appeared in more than two pointings: to accumulate the maximum possible 
number of two-point matches, we match each appearance of a source with every other appearance. 




tion, which is dictated by the radiometer equation: where Tgys is the system temperature, Na is the num- 

ber of antennas, r is the integration time, and Av is the 
bandwidth. The term {2kTsys / Aeff) is also known as 
AS = , , (10) the SEFD. (Note that the actual source count distribu- 

te// V^M^Vl^^l)^ 
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Fig. 7. — The primary beam's FWHM 0i/2 vs. observing frequency, shown over the relevant portion of the ATA's 0.5-10 GHz frequency 
range. Points on the solid line are the expected values from Equation [S] and the two large points are the FWHM values that we calculated 
using the method, as reported in i|4.2l 



tion as a function of flux is more complicated, and is 
only proportional to 5"^ over the ~1-1000 mJy range of 
fluxes detected in the PiGSS survey.) 

The RMS noise should decrease with an increasing 
number of antennas. This fact, combined with the in- 
creasing number sources that will be detectable as RMS 
decreases, should result in better estimations of the 
beam's FWHM using the fitting method we describe 



in l2.2l To test our prediction, we simulate datasets from 
telescopes with an increasing number of dishes, to which 
we then apply our fitting method in the same way that 
we fit the real ATATS and PiGSS data. We assume that 
the arrays in our simulations are not confusion limited. 

We simulate observations from seven different arrays 
with numbers of antennas ranging from 42 to 2688 in 
powers of two. For each of the seven arrays, we generate 
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Fig. 8. — The percentage uncertainty in the FWHM of the tele- 
scope's primary beam as a function of antenna number. The solid 
line is a power law fit to the simulation results, and has an index 
of 1. 



and analyze 1,000 different datasets. 

For each simulated dataset, we generate a field of the 
appropriate number of sources by distributing them ran- 
domly throughout a 12.6 deg^ field identical to the seven- 
pointing PiGSS field shown in Fig. [2j We then add Gaus- 
sian noise with a mean equal to the RMS of each source's 
flux density. After distributing the sources across the 
field, we "observe" all of the sources in each of the seven 
pointings, where the observed flux density of the source 
is equal to the intrinsic flux density scaled down by the 
gain as a function of distance from the pointing center. 
As we have throughout this work, we assume the pri- 
mary beam is a circular Gaussian with a FWHM equal 
to 1.10°, which we calculated in after analyzing the 
PiGSS data. We eliminate any observed sources whose 
fluxes have been scaled down to less than the RMS, and 
we only count 5a sources as detections. After accumu- 
lating the lists of sources in each of the seven pointings, 
we use the same procedure that we use in §2.21 to deter- 
mine the best-flt FWHM of the primary beam and its 
uncertainty. 

We flnd that as the number of dishes increases, the 
uncertainty in the primary beam's FWHM decreases like 
a power law with an index of ^^1 (see Fig. [8|). The first 
data point, corresponding the ATA-42, reports a FWHM 
uncertainty of 0.03°, which is consistent to within a factor 
of three with the results from the PiGSS data reported in 
i j4.2l As expected, our results return a median FWHM 
equal to 1.10°, which is equal to the value we used when 
we simulated the data; and each simulation returns a 
xi value near 1, which is reasonable since we haven't 
introduced any systematic errors into our calculations. 

The decreasing error in the calculated FWHM of the 
beam with increasing dish number will be important for 
future telescopes with large numbers of dishes . For ex- 
ample, with the SKA-3000 (Schflizzi et al."2007), one wiU 
be able to characterize the primary beam much better by 
using images taken with a fraction of the integration time 
characteristic of the above ATA-42 observations. Extrap- 
olating the power law in Figure [5] to an array with 3000 
antennas, we can see that in the absence of systematic 
errors one could estimate the FWHM of the SKA-3000 
to within 0.02%. 



6. CONCLUSIONS 

While consistent with our minimization results, the 
method of calculating the beam's FWHM using the two- 
point method has systematic problems, including the fact 
that the best flux density-distance pairs lead to the least 
reliable FWHM values. We also flnd the minimization 
to be superior to the two-point method because of its lack 
of systematic problems, its ability to use all of our data in 
the flt, and its compatibility with more complex models. 

Overall, both our analysis of real ATA data and our 
simu lation confirm th a t the theoretical FWHM values 
from I Welch fc DeBoed (l200l) are reasonable, as they are 
consistent with our results at the frequencies of both 
the ATATS and PiGSS surveys. The theoretical values 
are also consistent wi th the results from radio hologra- 
phy (iHarp et al.ll2010D and off-axis in easurement of point 
sources (jMacMaho n fc Wrightl [20091) . 

The results of our simulation, which calculates the un- 
certainty in the FWHM of the primary beams of tele- 
scopes with different numbers of dishes, will be impor- 
tant for the primary-beam characterization of future 
LNSD telescopes such as the ATA-350 and the SKA. 
The dramatically increased number of sources detectable 
by these telescopes will allow much better constraint of 
beam properties, thus improving the accuracy and fi- 
delity of mosaicked images. Cornwcll et al. (1993) note 
that in the absence of systematic errors, one can achieve 
a dynamic range in a mosaicked image of '^3,000 and a 
fidelity index of ^-^150 by knowing the primary beam to 
an RMS error of ^1%. We can measure the beam shape 
of the ATA to within a tenth of a percent, suggesting that 
we have the potential to achieve extremely high dynamic 
ranges and fidelity indices in ATA mosaics. 

An extension of our model could be developed to al- 
low fitting not only of the beam's FWHM, but also of 
its ellipticity and its angle on the sky. Additional effects 
such as spillover from the ground, mechanical deforma- 
tion of the dish, frequency and baseline dependence of the 
beam, and time variability of sources are surely factors 
to be considered; however, most of these are probably 
relatively minor effects when compared with the effect 
of using a simple, Gaussian model for a beam that we 
know is non-Gaussian at some level. Regardless, all of 
these effects should be taken into consideration in future 
primary beam studies on telescopes such as the ATA-350 
the SKA, and other radio interferometers. 
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